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Abstract 



We perform a detailed numerical study of energy-level and wavefunction 
statistics of a deformable quantum billiard focusing on properties relevant 
to semiconductor quantum dots. We consider the family of Robnik billiards 
generated by simple conformal maps of the unit disk; the shape of this family 
of billiards may be varied continuously at fixed area by tuning the parame- 
ters of the map. The classical dynamics of these billiards is well- understood 
and this allows us to study the quantum properties of subfamilies which span 
the transition from integrability to chaos as well as families at approximately 
constant degree of chaoticity (Kolmogorov entropy). In the regime of hard 
chaos we find that the statistical properties of interest are well-described by 
random-matrix theory and completely insensitive to the particular shape of 
the dot. However in the nearly-integrable regime non-universal behavior is 
found. Specifically, the level-width distribution is well-described by the pre- 
dicted distribution both in the presence and absence of magnetic flux when 
the system is fully chaotic; however it departs substantially from this behavior 
in the mixed regime. The chaotic behavior corroborates the previously pre- 
dicted behavior of the peak- height distribution for deformed quantum dots. 
We also investigate the energy-level correlation functions which are found to 
agree well with the behavior calculated for quasi-zero-dimensional disordered 
systems. 

PACS numbers: 73.40.Gk, 05.45.-hb, 72.20.My, 73.20.Dx 
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I. INTRODUCTION 



The relevance of concepts from the theory of "quantum chaos" to mesoscopic physics has 
become increasingly clear as nanostructure technology has achieved controlled fabrication of 
systems smaller than both the elastic and inelastic scattering lengthlll. Quantum chaos is the 
generally accepted term for properties of a quantum system associated with classical chaos 
(or the classical transition to chaos) . In the most recent high- mobility nanostructures trans- 
port is ballistic and the dominant scattering mechanism is the reflection of the electrons at 
the boundaries of the structure which (depending on the nature of the confining potential) 
may generate classically chaotic, mixed or integrable dynamics. Thus they present experi- 
mental possibilities for the application and testing of concepts from the theory of quantum 
chaos in condensed matter physics. Disordered mesoscopic systems also are undoubtedly 
chaotic classically, and recent work has emphasized the similarity between disordered quan- 
tum systems and the ballistic systems which have chaotic boundary scattering. In particular 
at low temperature both type of systems exhibit sample-specific mesoscopic fluctuations in 
various physical properties as a function of external parameters such as magnetic field. The 
ballistic systems differ from the disordered systems however in two ways. First, as we will 
use the term, a disordered system generates elastic scattering of electrons on a scale / which 
is short compared to typical sample dimensions L. This means that transport is diffusive 
on a scale smaller than the system size and because the diffusion process is dependent on 
dimensionality the statistical properties in general depend on the spatial dimension. Chaotic 
ballistic systems on the other hand have no relevant transport length smaller than the system 
size and thus many of their properties are insensitive to the spatial dimensionality (they are 
said to be "quasi-zero-dimensional" ) . It turns out that this difference leads to differences in 
the statisticalproperties of disordered and chaotic systems in certain regimes of energy and 
temperatureBcl. Second, there is reasonable evidence that ballistic systems may be fabricated 
with geometries and potential profiles which generate nearly-integrable classical dynamics, 
thus it becomes worthwhile to consider models which describe the transition to chaos and 
not just fully chaotic dynamics. 

There are two main approaches to the quantum theory of chaotic systems. The approach 
through semiclassical quantum mechanics pioneered by Gutzwilleri, and the approach based 
on the theory of random matrices first applied to quantum chaos by Bohigas, Giannoni and 
Schmiti. The former approach makes a more direct connection to the classical mechanics 
and has had major successes recently in atomic physics. However the confining potential in 
the microstructures studied experimentally is rarely well enough known to justify theoreti- 
cal work relying on specific classical orbits. Instead either a wholly statistical approach, or 
a combination of semiclassical and statistical ideas has been applied to quantum chaos in 
mesoscopic systems. Three measured physical effects which have been proposed as mani- 
festations of quantum chaos in mesoscopic transport are: 1) The resistance fluctuations in 
GaAs quantum wires coupled strongly to an electron cavityS'i^. 2) The weak localization 
effect in the same system&i. 3) The fluctuations in the Coulomb blockade conductance 
peaks in quantum dots weakly coupled to leads&H. The electron cavity conductance 
fluctuations have been described by a combination of statistical and semiclassical theory, and 
most recently certain properties have been derived from random-matrix theorylll0. Whereas 
the Coulomb blockade peak fluctuations in quantum dots have been described completely 



2 



statistically, using only random-matrix theory. The quantum dot conduction experiments 
are analogous to strongly resonant scattering in atoms or nuclei for which the properties of 
a single quasi-bound state can be probed. In micron-size semiconductor quantum dots at 
the typical experimental density it is estimated that the single-particle level-spacing (or the 
excitation energy to the first excited state) Ae ~ 0.05 meV ~ 500 mK and therefore these 
systems may be studied in the regime kT < Ae where indeed only a single quasi-bound 
state participates in the resonance. 

In this paper we will concentrate on the statistical approach to quantum chaos through 
random-matrix theory. This theory was originally developed by Wigner, Dyson and others to 
explain statistical properties of compound nuclear resonancesBEl but later it was conjectured 
by Bohigas-Giannoni-Schmiti to describe any quantum system whose classical analogue is 
fully chaotic. There now exists substantial numerical evidence supporting this conjecture 
as well as an analytic argument due to Berryll which applies to a particular statistical 
property (known as the A3 statistic) measuring the long-range rigidity of the spectrum. 
In RMT complex systems are represented by ensembles of Hamiltonians with statistically- 
independent matrix elements. These ensembles can be shown to have maximum statistical 
entropy (subject to a small number of constraints)ill. There exist three different symmetry 
classes of such ensembles characterized by the number (3 = 1, 2, 4 of independent components 
of the matrix elements in the Hamiltonian, H: In zero magnetic field without spin-dependent 
scattering H is real, (3 = 1, and the corresponding ensemble (invariant under orthogonal 
transformations) is denoted the Gaussian orthogonal ensemble (GOE). In non-zero magnetic 
field H is complex, (3 = 2, and the ensemble (invariant under unitary transformations) is 
denoted the Gaussian unitary ensemble (GUE). Finally for strong spin-orbit scattering H 
is quarternion real, (3 = 4, and the ensemble (invariant under symplectic transformations) 
is denoted the Gaussian symplectic ensemble (GSE). In this paper we shall neglect spin 
effects (since spin-orbit scattering is neglibile in the ballistic microstructures) and hence 
only consider applications of the GOE or GUE. 

One well-known property of these statistical ensembles is the appearance of strong level- 
repulsion due to the lack of any conserved quantum numbers (other than the energy). This 
makes the probability density of level spacings, s, tend to zero as s — > unlike the Poisson 
distribution of uncorrelated random variables (which has been shown to describe typical 
integrable systemsii). Using RMT it is possible to calculate the distribution Pf3{s) in the 
spectrum normalized such that the local average level spacing is one. For RMT of two- 
dimensional matrices the level spacing distributions are 



Pi{s)= fs exp(-fs2) (GOE) (1) 
P^is) = f .2 (_i,2j (GUE) (2) 

Note that this distribution tends to zero as as s ^ 0. Although not exact for N x N 
random matrices, these formulae (known as the Wigner surmise) are excellent approxima- 
tions to the exact results and are conventionally used for comparison to statistical data. The 
appearance of the RMT spacing distributions at the classical transition to hard chaos has 
been verified in many numerical studies (although an analytic derivation has still not been 
found). Hence in this work we will only use the spacing distribution as a diagnostic for the 
applicability of RMT. 
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The theory of these random-matrix ensembles has tended to focus on spectral statistics 
and not on statistical properties of the eigenstates; however in the application to quantum 
dot resonances treated below we shall find that it is the eigenstate statistics which are most 
easily measured. Very recently a great deal of progress has been made in the study of 
ensembles in which the random-matrix varies as a function of an external parameter such as 
magnetic flux. Various correlation functions of the spectra have been calculated and shown 
to be universal upon rescalingil'B. These correlation functions are not easily measured 
in current experiments on quantum dots, but may be accessible with some effort!!. Since 
such correlation functions are only calculated analytically for random ensembles it is again 
worthwhile to test their applicability to a given chaotic dynamical system. We perform 
detailed comparisons of this type below. 

Although the applicability of random-matrix theory to quantum systems which are clas- 
sically chaotic is now reasonably well established, most dynamical systems are neither fully 
chaotic nor fully integrable, but instead have a mixed classical phase space described topo- 
logically by the Kolmogorov-Arnold-Moser (KAM) theoremS and related results. The ex- 
pectation is that RMT will give valid and universal results for fully chaotic hamiltonians 
and will break down in some manner as regions of stability appear in the phase-space. It 
is of interest then to look at a model for an ensemble of quantum dots which can span the 
range from integrability to hard chaos, and can also test universality in the chaotic limit. A 
model of this type was introduced by RobnikHil and generalized by Robnik and BerryH. 
In this model the quantum dot is represented by a deformed circular well with infinite walls; 
the deformation is described by a quadratic or cubic conformal transformation. This model 
has a number of attractive features for our Durposes. 1) The classical mechanics is well- 
understood and has been studied in detaili3§. It has been shown that as the parameters 
of the conformal transformation are varied the well can undergo a standard KAM transi- 
tion to chaos. In addition we show that if the parameters are varied in a different manner 
a sequence of different chaotic billiards are generated with roughly equivalent degrees of 
classical chaos. Thus one can use this set for statistical averaging. 2) Robnik showed that 
a very efficient numerical algorithm exists for obtaining a large number of eigenstates and 
eigenvalues for this model. 3) Berry and Robnik showed that an Aharonov-Bohm flux may 
be simply introduced into the well to break time-reversal symmetry (without changing the 
classical dynamics) and that this would only introduce minor changes in the numerical al- 
gorithm for solution of the Schrodinger equation. Thus the model is suitable for testing the 
universality of statistical properties in the chaotic regime and their possible breakdown in 
the mixed regime for both the orthogonal and unitary ensembles. Preliminary results of this 
investigation of significantly less breadth have been published elsewhere^. 

The paper is organized as follows. In Sec. II the deformable billiard is defined in terms 
of a conformal map. In Sec. Ill the classical dynamics of the billiard is discussed with the 
emphasis on characterizing the degree of chaos as a function of deformation. In Sec. IV the 
quantum mechanics of the billiard is discussed, and the single-electron Schrodinger equa- 
tion is solved with three different boundary conditions (Dirichlet, Neumann, and general), 
and the respective level statistics are calculated. In Sec. V using the shape as the external 
perturbation the energy level correlation functions are calculated for the normalized spec- 
tra and comparisons are made with conjectures for their form stated in the literature. In 
Sec. VI we turn to the problem of fluctuating Coulomb blockade conductance peak heights 
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in quantum dots. A model of a Coulomb blockade device is made based on the deformable 
billiard, and the peak height distributions are calculated and compared to those obtained 
from random-matrix theory. In Sec. VII some experimental implications are discussed and 
a new experiment is suggested aimed at testing our results. Finally in Sec. VIII we present 
a concluding discussion. 



II. DEFINITION OF THE DEFORMABLE BILLIARD 

A mathematically simple way of defining a continuous deformable family of billiards was 
introduced by RobnikS and Berry and Robnikil. It is based on a conformal mapping of the 
unit disk. To make this paper self-contained and to facilitate the presentation of some new 
developments of the method we shall briefly describe the technique below. 

As shown in Fig. |I] we study the simply connected domain T) with a possibly irregular 
shaped boundary dV in the uv plane. The open interior of the domain is denoted V\&D. 
The deformed billiard is defined by the infinitely hard wall potential V{u,v) satisfying 

V(u.v) = {° ■ <"'"!^^\^ (3) 

y oo , {u,v) ^ V. ^ ' 

Following Robnik0 we define the shape of P by a conformal mapping w of the unit disk C 
in the xy plane to T> in the uv plane. Parts of the treatment are more conveniently carried 
out by introducing the complex coordinates z and w. 

z = X + iy w = w{z) = u{z) + iv{z). (4) 

The boundary &D in the w plane is thus given as the image w{C) of the boundary of the 
unit disk dC in the 2;^lane. Throughout this work we study the cubic mapping introduced 
by Berry and Robnil< 



J29 



, , z + bz'^ + ce^^z^ ^ , , 

w(z) = , z eC, (5) 

where b, c, and 6 are real parameters chosen such that > for ^ G C Two sequences 

of deformed billiards are show in Fig. |^. This cubic form of w{z) is the simplest conformal 
map resulting in a billiard with no spatial symmetries (see Fig. ^d). The even simpler 
quadratic map (c = 0) generates a family of billiards with reflection symmetry (see Fig. ^); 
Robnik and Berry pointed out that such spatial symmetries can prevent a magnetic flux from 
generating the orthogonal to unitary transitioned hence we maintain the more general form 
although the simple quadratic case is sufficient for many of our calculations. The deformed 
billiard is given by: 

dV = {wiz) : 1^1 = 1} 
V\dV = {w{z) : \z\ < 1} 

The normalization in Eq. (^) insures that the area of V remains vr for any value of the 
parameters. 
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For later reference we give the explicit form of the Jacobian in the polar coordi- 

nates (r, 6) of the xy plane: 

\w'{z = re^^) p = [1 + + 9cV^ + 6cr^ cos(20 + 6) 

+46r cos(6') + 12bcr^ cos{9 + 5)] 
/(l + 262 + 3c2). (7) 



III. THE CLASSICAL DYNAMICS OF THE DEFORMED BILLIARD 



Robnik conducted a thorough studycZl of the classical mechanics of a point mass moving 
freely in the deformed billiard for the quadratic case (c = 0, see Fig. |^). More recentlyE^ 
Hayli et al. have extended his results. Thus in contrast to many previous works on bal- 
listic microstructures in which a discretized version of the Schrodinger equation has been 
studiedB'iS in this case one has a detailed knowledge of the classical dynamics of the rele- 
vant quantum system. Robnik showed how as a function of b, starting from 6 = 0, the system 
evolves according to the KAM theorem into a mixed phase-space exhibiting soft chaos and 
eventually to fully developed chaos. The main tool for determining the degree of chaos has 
been to construct the Poincare surfaces of section for the bounce map and to calculate the 
Kolmogorov entropy (we will define these below). The last large islands of stability in the 
Poincare section disappear around b = 0.25 (the value at which the billiard ceases to be 
convex) and Robnik originally conjectured that the transition to hard chaos occured at this 
valueEH. However recent work has shown that very small islands of stability spawned by the 
bifurcation of the final stable two-cycle persist up to 6 ~ 0.28, but the precise value of b 
at which hard chaos sets in is not known. However it has recently been provenil that the 
quadratic billiard at 6 = 0.50 is fully chaotic making this only the third billiard (along with 
the stadium and Sinai billiards) for which hard chaos has been demonstrated analytically. 
In practice the islands of stability have negligible weight for b > 0.25 and the statistical 
properties of both the classical and quantum mechanics are consistent with fully developed 
chaos. 

The Poincare section is constructed by plotting the phase space coordinates of the particle 
each time it reflects specularly from the boundary of the billiard. The Poincare section 
is area-preserving if it is constructed using conjugate canonical variables; a useful set of 
coordinates for billiard^ is the arc length a, measured from the origin to the point the 
particle hits the boundary, and the tangential momentum sin(x) at this point (x denotes the 
angle of incidence). In Fig. ^ we show Poincare sections corresponding to the deformation 
sequence of Fig. |^. For simplicity here we replace a, the arc-length, by 6 the angular 
coordinate. For small deformations these are approximately proportional and since cr is a 
monotonically increasing function of 6 the topology of the surface of section in 6 is the same 
as in a. (strictly speaking our use of 6 introduces some small violation of the area-preserving 
property which is unimportant since we are only using the SOS for illustrative purposes). 
The surfaces of section in Fig. ^ clearly show how the deformed billiard gets more chaotic 
as the deformation is increased from zero. The many large stable islands at the beginning 
of the sequence shrink and vanish sequentially until one obtains the featureless surface of 
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section plot shown in Fig. In contrast, if one starts from a deformed billiard as in Fig. ^3 
and progresses through further deformations, all the surface of section plots look like the 
featureless plot in Fig. indicating that all of these billiards are close to hard chaos. 

Quantitative information on the degree of chaos can be obtained by calculating the 
Lyapunov exponent and the related Kolmogorov entropy00. Following Benettinii the 
Lyapunov exponent can be estimated numerically by examining two trajectories starting 
at the phase space points = {ctq, sin(xo)} and = {ctq, sin(xQ)} separated by a small 
distance do = — ~ 10^^°. After one bounce the two points $,1 and ^\ are reached, 
the distance di between them is calculated, and the reduction ratio /3i = d^/di is formed. 
A new rescaled point ^5* is defined by = + A (^5 " ^i)- Starting from and ^5* a 
second bounce is calculated, and from ^2 and ^2 ^ new reduction ratio /?2 and a new rescaled 
point ^2* is formed and the process is repeated. The Lyapunov exponent A which generally 
depends on the initial phase space point can then be calculated as 

AK,sin(xg)) = - jim (^EI^a) , (8) 

where t„ is the traversal time or the total length in real space of the total n-bounce orbit. 
Prom the local measure of chaoticity constituted by A, which depends on the initial points 
in phase space, a global measure the Kolmogorov entropy h can be formed: 

where C is the perimeter of the billiard and (cr, sin(x)) is the initial phase space point. The 
Kolmogorov entropy h provides a quantitative measure of the degree of chaos in the billiard. 
In Fig. I /i is shown for the deformation sequences of Fig. 0. Note how h grows monotonically 
for the sequence (a), while it is fairly constant and significantly bigger for the sequence (b). 
Thus the Kolmogorov entropy calculations indicate that the deformation sequence Fig. ^3 
is "more chaotic" than the quadratic billiard even at 6 = 0.50 where it is known to be fully 
chaotic. Of course the K-entropy is essentially an average property of the phase-space and 
does not allow us to exclude some very small regions of stability in the sequence Fig. ^p. 



IV. QUANTUM MECHANICS OF THE DEFORMED BILLIARD 

Following Berry and Robnikil one can calculate the quantum states of a single particle 
moving in the interior of the deformed billiard given by Eq. (^). The billiard is threaded by 
an Aharonov-Bohm flux tube of strength a$o through the origin of the uv plane. Here $0 
is the flux quantum h/e and a is a dimensionless real number. Choosing the gauge 

the Schrodinger equation in w coordinates, 

^(-zW^ + eA{w)fm{w) = E^H, (11) 
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becomes the following in the polar coordinates (r, 6) of the xy plane: 

-f^^(r,^) +e|u7'(re^^)|2^(r,^) = 0, ^ 

where now the energy e is given in units of fi^ /2mB? and lengths in units of i?, R being 
the radius of the circle C. The spectrum of the Hamiltonian is periodic in a with period 
1. Except for the two values a = 0, 1/2 for which a real nondiagonal representation of the 
Hamiltonian can be found the Hamiltonian is a full complex Hermitian operator. Hence 
GUE statistics is expected for all a except for values close to and 1/2 at which GOE 
statistics is anticipated. Below whenever we wish to focus on the GUE case we choose the 
value a = 1/4 to stay as far away as possible from the GOE values. 

Thus using the fact that the deformed billiard is obtained by a conformal map of the unit 
disk we may replace the original problem on the irregular domain T) by an equivalent problem 
on the unit disk moving under a rather simple "potential" proportional to |w'(re*^)p. In 
the simplest cases of Dirichlet or Neumann boundary conditions discussed in sections A 
and B below the boundary conditions of the equivalent problem are identical to those on 
the original domain. Since a convenient basis which satisfies these boundary conditions is 
available (products of Bessel functions, Jui'jr), and exponentials, exp{iW)) Eq. ([12|) may be 
expressed as a linear equation in that (infinite) basis and after truncation may be solved 
numerically. Typically in this study we truncate Eq. (0) to the subspace of the 1000 lowest 
energy Bessel solutions. As noted by Robnik, the dependence of this linear equation on 
the shape parameters b, c, and 6 can be factored out allowing a very efficient numerical 
method of solution as the shape is varied. To discuss physical properties relating to isolated 
quantum dots it is reasonable to consider this model with Dirichlet boundary conditions; 
however we are also interested in fluctuation properties related to conductance. In this case 
the deformed dot must be coupled to leads (typically through tunnel barriers) and more 
general boundary conditions must be imposed to describe its scattering resonances. Hence 
in the following subsections we discuss the solution of Eq. (p!^ ) under Dirichlet, Neumann 
and flnally general inhomogeneous boundary conditions. 



A. Dirichlet boundary conditions (DBC) 



Eq. (12) with the Dirichlet boundary conditions (DBC) 

*(w) = 0, for w e dV. (13) 



was considered in the original work of Berry and Robnikc^. As just noted, a simplifying 
property of the conformal mapping is that if \l/(w) satisfy DBC on the domain D, so will 
'^{w{z)) on the circle C. We therefore expand the solutions of the Schrodinger equation (|12|) 
with DBC in terms of the eigenstates |0i,n) of the analogous Schrodinger equation for the 
unit disk C also with DBC. These are given by 

{r,0\<Pl,n) = ^j, ^, rJ|i_,|(7K-a|,nr)e''', (14) 
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where r < 1, J,y is the Bessel function of order z/, 7i,^„ is the nth root of J^, and / is an 
integer. The eigenvalues of |0i,n) are given by 

-VV^n)=7'-.|l0^n). (15) 

Note that degeneracies occur in the energy spectrum in case of zero AB-flux (a = 0). 

It is convenient to enumerate the basis states not according to the quantum numbers 
(/, n) but according to the position j of the eigenvalues 7^ arranged in ascending order. 
The eigenstates ipp of the Schrodinger equation in the deformed billiard with DEC are now 
written as the superposition 

\i^p)=Npj:^\<p^)^ (16) 

j=l 1] 

where Np is a normalization constant. The particular form of the expansion coefficients is 
chosen so that the transformed Schrodinger equation (|12D becomes a Hermitian eigenvalue 
problem. The coefficients c'-^-' and the eigenvalues Sp are found from Eq. ([1^) by insertion 
of the superposition Eq. ([T6|) followed by multiplication with (r, 6'|0j) and integration over 
r and 9. After rearrangement the eigenvalue problem then takes the form 

M.,cf = icf , (17) 
The explicit form of the Hermitian matrix M using Eq. (0) for the Jacobian is the following: 



+ 5i- 

+ 5l. 

' ''7. 1 



,,MI^? + ^<^I^) (18) 
,,,6ce+^^/(;V(l + 26' + 3c2), 



where the integrals llf' are given by 

7.7,^;(7.)^;(7,) 



(19) 



with z/j = I /j — a I . For given values of the parameters b, c, and 6 we construct the truncated 
M-matrix for the 1000 lowest eigenstates \4>j)- The 300 highest eigenvectors of M corre- 
sponding to the 300 lowest eigenstates of our original Hamiltonian are accurate enough for 
our analysis. Typically the eigenvalues Sp are more accurate than 0.001 times the average 
level spacing, and among the 300 levels only a few of the highest ones, where the presence 
of the matrix truncation is felt most strongly, have an accuracy as low as around 0.1 times 
the average level spacing. 

Once the energy spectrum is obtained for the lowest 300 eigenstates various level statistics 
can be calculated and compared with the predictions of RMT. For specific shapes (parameter 
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values) this has been done previously by Berry and Robnikc3 and later by Goldberg et alM, 
finding excellent agreement with e.g. the spacing distribution predicted by RMT in the 
chaotic regime. Since we will use a range of shapes not covered by these studies we confirm 
the expectation that these results extend to the wider class of chaotic billiards studied below. 

The first step in the energy level analysis is to "unfold" the spectrum. This procedure 
corrects for the average variation in the density of states with energy so that the mean level- 
spacing of the unfolded spectrum is always unity and data in different energy intervals may 
be merged. For a smooth billiard obeying DBG the Weyl formula states that the average 
number of states with energy less than e, N{e), is given (for high energies) byS: 

where C is the perimeter of the billiard, which we calculate by numerical integration. The 
unfolded eigenenergy corresponding to the eigenenergy ej is given by N{ej), and the jth 
unfolded level spacing sj is thus 

s^^Nie,^,)-NiE,). (21) 

After unfolding the calculated energy spectra a histogram for the distribution of the level 
spacings can be constructed, an example of which is shown in Fig. |^. The hypothesis that 
the histograms are following the GOE distribution of Eq. (|l]) or the GUE distribution of 
Eq. (0) can now be tested statistically. The quantitative measure we use is the one-sided 
X^-test, which estimates the probability -P(x^) for finding a x^-value bigger than the actual 
observed one. 

In Fig. ^ we show the results of the x^-test for the GOE and GUE distributions. We have 
calculated the energy spectra for 78 different values of 6 in the interval [0, vr] for constant b 
and c, both being equal to 0.2 (the deformation sequence shown in Fig. Qd). The billiard 
was threaded by an AB-flux tube of strength | of a flux quantum. It is seen, as expected 
from Ref. that the level distribution with statistical significance is given by the GUE 
distribution in most of the parameter range. However, we also note that for 5 ~ and 
(5 PS TT where the billiard is almost symmetrical around the w-axis, the GOE distribution 
better describes the spectrum. When comparing the results for the classical mechanics 
in Sec. Ill with those obtained here for the quantum mechanics, we obtain yet another 
verification of the observation that the energy level distribution of a quantum system whose 
classical counterpart exhibits hard chaos follows Wigner's distribution. Although this result 
was expected, we note that our study provides a more complete test of the validity of this 
statement than previous studies in that we test a range of billiards with roughly constant 
Kolmogorov entropy. 



B. Neumann boundary conditions (NBC) 

A central issue in this paper is to study not completely isolated but nearly isolated 
systems where the electrons can escape. Therefore we have to abandon the DBG of Eq. (|T3|) 
which are relevant only for closed systems. As we will discuss later, a partially open system 
will require not just boundary conditions but also matching conditions for solutions inside 



10 



and outside. Such matching conditions can be expressed in terms of any basis set for 
the region V which does not cause the wavefunction to vanish identically (as do DBC). 
The mathematically simplest alternative choice is Neumann boundary conditions (NBC) for 
which the normal derivative vanishes everywhere on the boundary: 

n-W(u7) = 0, w G dV, (22) 

where n is the outward pointing normal of the boundary at the point w. 

Although not treated in previous work NBC maintain the simplicity of DBC in that 
they are preserved by the conformal map, i.e. if '^{w) obeys Eq. (^) so will '^{w{z)) for 
z e dC, and there exists a simply basis for the unit disk satisfying NBC. We can therefore 
solve the Schrodinger equation Eq. (0) with NBC following the steps in Sec. IV-A. The 
only modification is a change of the radial part of the basis functions defined in C as they 
appear in Eq. ( P^ ) to the following (where a tilde is used to distinguish NBC from DBC): 



{r. e\4>,) = ,^ /^'/^J^ — J|,^_.|(7,r)e^'^-^ (23) 

We have again ordered the basis set in ascending order according to the eigenvalues 7^ , with 

7j now being a root of the derivative J'^. of the Bessel function Ju^. The solution of 
Schrodinger equation with NBC in the deformed billiard is written in analogy with Eq. (|TBp 
as: 

_ 00 ~ip) 

\^p) = NpJ2^\^,), (24) 
j=i 7j 



M^r^f = Igf , (25) 



The matrix equation ([T7| ) becomes 



The explicit form of the Hermitian matrix M is obtained from Eq. (|l^) by substituting the 



integrals I-'^^ with the integrals I^j'^ given by 



(26) 



where z/, = |/j — a|. The energy level statistics under NBC can now be studied by calculating 
the normalized energy spectrum using numerical diagonalization of the truncated 1000 by 
1000 sub-matrix of M. The Weyl formula for the NBC case isH 

where in comparison with the DBC case Eq. ( PPD the only change is the sign of the second 
term on the right-hand side. 
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Although the change from DBC to NBC leads to relatively minor changes in the method 
of solution, the relation between the quantum mechanics and the classical mechanics studied 
in Sec. Ill becomes less clear since it is natural to associate the hard walls in the classical 
problem with DBC in the quantum problem. Thus it is of some interest to see if the relation 
between classical chaos and RMT is independent of this change in the boundary conditions 
of the quantum problem. The histograms of the level spacing distributions for the NBC 
case are calculated for the same AB-flux and the same shape deformation sequence as for 
the DBC case, and x^-tests for the GOE and CUE level statistics are performed. The 
result is shown in Fig. |^, where it is seen that again the level spacing distributions are well 
described by CUE statistics, except near the special values of a = 0, 1/2. We observe that 
the crossover region between GOE and CUE statistics near the 5 = is wider in the NBC 
case than the DBC case. We have no physical explanation for this behavior and suspect 
it is a numerical artifact. Nonetheless the results overall indicate that the correspondence 
between classical chaos and RMT behavior for the spacing distribution is independent of 
the boundary conditions. 



C. General boundary conditions (GBC) 

As noted in our introduction, one proposed application of RMT to quantum dots relates 
to the peak amplitude fluctuations in resonant tunneling through the dot. In such a case 
electrons may enter and leave the dot in certain directions through leads separated from 
the dot by tunnel barriers. Thus the problem of interest has the geometry illustrated in 
Fig. 1^. The Schrodinger equation now has solutions at all energies although resonances 
will still occur at an energy-spacing comparable to that of the closed system. For such a 
geometry the physically relevant boundary conditions are Dirichlet on the portion of the 
billiard unconnected to the leads (we denote this part of the boundary dVi) combined with 
matching conditions at the leads. In general one may expand the wavefunction inside the dot 
in any basis which allows the matching procedure at the leads, (i.e. such that the sum over 
basis functions is uniformly convergent). Unfortunately Dirichlet boundary conditions at 
the leads do not satisfy this condition, thus one is forced to use a basis which satisfies mixed 
boundary conditions on dT>] DBC on &Di and general inhomogenous boundary conditions 
on the complement dT>2. 

5f(w) n- W(u7) + ^(u;) = 0, w e dV 

g{w) = 0, we dVi (28) 
g{w) 7^ 0, we dV2 

We refer to these as general boundary conditions, GBC. Unlike DBC and NBC, GBC 
are changed by the conformal mapping; in fact the function g{w) will transform to 
g{w{z))/\w'{z)\. However a more fundamental problem is that no simple basis set exists 
on the unit disk which satisfies GBC, so no straightforward extension of the DBC and NBC 
approach is possible. Of course the problem we are encountering occurs in other contexts 
(e.g. electromagnetic waveguides) and with significantly more effort we are able to use ex- 
isting techniques to get an approximate solution. We employ a two-step procedure. First we 
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find the eigenstates of tlie deformed billiard obeying DBC, and then we use the method de- 
veloped by FeshbachS to perturb the DBC at &D2 with the GBC of Eq. (|2^) . The strength 
of the perturbation is given by the function g{w). If g{w) were zero everywhere along V we 
recover DBC, and if g{w) were infinite everywhere along V we get NBC. The method we 
use below only applies in the case of zero magnetic field. 

Our point of departure is the eigenstates {ipp) of Eq. (0) obeying DBC. We have 

V%iw) + Spi/jpiw) = ijjpiS) = 0, (29) 

where 5* here and in the following denotes a point on the boundary &D. Our aim is to find 
a solution Xx in V satisfying GBC: 

V^Xx{w) + ExXx{w) = 
giS)n-VX,iS)+X,iS) = 0. (30) 

We begin with the Green's function Gx{w\wq) which satisfies DBC: 

VlGx{w\wo) + ExGx{w\wo) = -S{w - wq) 

Gx{w\S) = 0. (31) 

Eqs. ([30| ) and (|3TD together with the use of Green's theorem yields 
Xx{w) = -fdS n-VoGxHS) Xx{S) 

JdV 

= + (f dSn-VoG{w\S)g{S)n-VXx{S). (32) 

JdV 

Taking the gradient Vu, of Eq. (|5^) and expanding VX^ in terms of V0p, 

VXa = E— (33) 

one obtains after some analysis^ the following eigenvalue problem involving the expansion 
coefficients SJ"^ and the eigenvalues Ex- 



E 



+ dSn-Vr,{S)g{S)n-V^,{S)\-^6,}4'^ = 0. (34) 

6p EpEq JdV ^ \ Ex } ^ 



Since we are expanding the wave function Xx obeying GBC in a basis set {ipp} obeying 
DBC some discussion of convergence properties is necessary. The expansion in Eq. (|33D 
converges only point-wise, and the convergence is slowest near the non-DBC part 81)2 of the 
boundary. There one might expect the occurence of "overshooting" (the Gibbs phenomena) 
of the summed series. However, the convergence is good in a least-squares sense, meaning 
that quantities involving integration over the whole domain T) converge well if the period 
of the smallest oscillation in the series is much smaller than the extension of the region 
where g{S) is non-zero. For semiclassical concepts to be relevant we require that the leads 
be many wavelengths across, but we also wish them to be small compared to the radius of 



13 



the dot. We chose to use leads with a width of about 7r/6 corresponding to about five times 
the smallest azimuthal wave length when the basis is truncated between 700 and 1000. The 
energies Ex, which is essentially the average of V^X;^ over V is an example of a quantity 
with a good least-square convergence. This convergence improves the smaller the strength 
for g{S). The DBC wavefunctions ipp appearing in Eq. ( |33D are only known in terms of the 
Bessel solutions of Eq. ([T6|) . It is thus convenient to map the line integral in Eq. ( ^If ) back to 
the unit disk where it becomes a sum of integrals along the perimeter of the unit disk with 
ipp and ipq replaced by Bessel functions. This mapping introduces an additional factor in 
the integrand A neat simplification can then be achieved by choosing the function 

g{S) as 

.c^^ _ . _ / 0.1 ,0e [0.5,1.0] U [4.5,5.0] . . 

g[b)-g[e ,9^ [0.5, 1.0] U [4.5, 5.0] ^'^^^ 

This choice cancels the factor mentioned above and allows simple analytical 

evaluation of the line integrals. This greatly improves the numerical tractability of the 
calculation. The strength of the perturbation g{S) is chosen so that the perturbation is 
strong enough to shift the individual levels of the order one mean level spacing and such 
that it is weak enough to avoid substantial mixing of the DBC-eigenstates invalidating the 
truncation we make at the second step of the procedure. In the following we show the result 
of calculating the energy spectrum in the GBC case based on the 700 lowest eigenstates of 
the DBC case. 

Once we have obtained the spectrum unfolding it requires a knowledge of the relevant 
Weyl formula since the subleading corrections depend on the choice of boundary conditions. 
Unfortunately there does not exist an analytical expression analogous to Eqs. (pOD and 
for GBC. However, the general theoryS implies that N(e) must have the form 

N{e) = ^e + By^ + C. (36) 

The constants B and C relating to the boundary and the curvature of the billiard can then 
be estimated by numerical fitting for each spectrum in question. We tested this procedure 
on the DBC and NBC spectra and obtained the known coefficients with only a few percent 
error. 

In Fig. P we show a comparison of tests of the numerical spacing distribution cal- 
culated with GBC compared to the Wigner surmise for the GOE and CUE for a range of 
spectra in the deformation sequence shown in Fig. for zero AB fiux. We see that GOE 
statistics is well-supported by the data. Taken together the results of sections A,B and C 
strongly support the conjecture that the correspondence between RMT spectral statistics 
and classical chaos holds independent of the boundary conditions on the quantum problem. 



V. PARAMETRIC ENERGY LEVEL CORRELATIONS 

Recently a great deal of theoretical work on disordered systems and RMT has examined 
the statistical correlation of energy levels when the hamiltonian varies as a function of some 
parameter such as fiuxiiH. The analytic results obtained apply to disordered systems 
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or to RMT ensembles not generated from a microscopic hamiltonian. Some nice numerical 
confirmation of the results has been obtained for ordered but irregular systems by Szafer 
and Altshuleril but not for models in which the classical mechanics was known. For example 
no numerical or analytic results are available for such correlations in the mixed regime. In 
this section we examine such correlations using the deformation parameters b and 6 as the 
control variable. We use the closed billiard with the simple DBC discussed in Sec. IV-A, 
and the deformation sequences shown in Fig. ^. In Fig. ^a only b is changed while in Fig. ^ 
only 6 varies, hence b and 6 will be the control variables relevant to each sequence, which 
we will denote generically as X. For a given sequence of parameter values the normalized 
energy spectrum is calculated. The resulting dimensionless energy levels are denoted ei{X). 
The prediction is that for fully chaotic systems certain correlation functions of the ei{X) 



are universal upon rescaling of X. Following Ref. ^we define the generalized conductance 
C(0) and the rescaled parameter 

'deiixy 



^^^^^\[~ix^ ^ ^ = VC(0)X, (37) 

where (■ ■ ■) is a statistical average over a suitable range of energy and/or X. The normalized 
energy level correlation function c{x) defined as: 

c(.) . (^^^) (38) 

is predicted to be a universal function which differs for GOE and GUE. 

The correlation function is calculated for the deformation sequence shown in Fig. ^3 
with the parameter 6 restricted to the interval [0.18,2.98] to avoid the symmetry points at 
and TT. The billiard is threaded by one quarter of a flux quantum and DBC is employed. 
The resulting normalized spectrum for the 50 levels between level number 225 and 274 



is shown in Fig. |10|. Since the spectrum changes both as a function of energy and as a 
function of S the spectrum has been divided into 5 by 8 boxes each containing 25 energy 
levels at 14 different values of S. Within each box the spectrum is fairly homogeneous, and 
the averaging of Eqs. (|37D and (^) is performed within each box. For each box 14 values 
of the pair {x,c{x)} can thus be calculated - the value of {0,c(0)}, however, is trivially 
{0, 1}. The resulting 560 points {x,c{x)} are then arranged in ascending order after their 
first component, they are grouped in 56 groups of ten, and the average value of each group 
is calculated. The final points are displayed in Fig. [ll|a. The same procedure is repeated for 
the case with zero AB-fiux and the resulting points are shown in the same figure. 



Fig. |Tl|a shows an excellent agreement between the predicted curves calculated in Refs. ^ 
and ^ for both the GOE and GUE cases. Thus our results support the conjecture that these 
correlations are universal and occur when the quantum system is classically chaotic. We can 
go further however and test whether these correlations actually coincide with classical chaos 
by looking at c{x) in the mixed regime. Our results shown in Fig. pT] b provide support for 
this stronger statement. c{x) in the mixed regime deviates strongly from its behavior for 
hard chaos. In particular a large dip appears in the correlation function which can be traced 
to the average distance to the first level anti-crossing. Unlike the chaotic regime where anti- 
crossings have large gaps, very small gap anti-crossings occur in the mixed regime leading 
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to a large parametric derivative of the energy levels. In fact it is knowrH that at such a 
anti-crossing in the mixed regime two wavefunctions "exchange identities" with the higher 
energy wavefunction having a spatial density close the the lower-energy wavefunction which 
has just been repelled. Thus we find that the spectral correlations are non-universal and 
inconsistent with RMT in the mixed regime. 



VI. STATISTICS OF COULOMB BLOCKADE CONDUCTION PEAKS 

As noted above, fluctuations in the spectrum of quantum dots are not yet easily accessible 
experimentally. The most striking fluctuation effect evident in the experimental data0'00 
are fluctuations in the peak height of Coulomb blockade resonances. These fluctuations 
reflect properties of the quasi-bound states(level- width fluctuations) and not of the spectrum. 
This contrasts with nuclear scattering resonances in which both snectral and level-width 
fluctuations are equally accessible. The reason for this difrerencdi& is that the quantum 
dot resonances correspond to the ground state energy of the system with N,N + 1, iV -|- 
2, . . . electrons and thus include the charging energy e^/C associated with the addition of a 
particle. Since this charging energy is approximately constant and is typically an order of 
magnitude larger than the single-particle excitation energy, Ae (or more precisely the energy 
to the first excited state for fixed N) the observed resonances are equally-spaced to a good 
approximation. In addition typically kT >> F (the mean level- width at zero temperature) so 
the resonances are thermally-broadend to a width ~ kT and only their amplitude reflects the 
level-width fluctuations. The amplitude fluctuations become maximal when kT < Ae and 
only the ground-state contributes to the resonance. In recent experiments00 Ae ~ 0.5 K so 
that this single-level regime is accessible. The cross-over between multiple-level and single- 
level tunneling leads to unusual and fluctuating temperature-dependences for the resonances 
until kT <C Ae, as was first understood by Meir et ahil. In earlier workEHl^ we have 
developed a detailed theory of the amplitude fluctuations in the single-level regime assuming 
that RMT describes the quasi-bound eigenstate fluctuations. Numerical tests of the theory 
agreed well for the GOE case but not as well for the GUE cas and were performed for a 
disordered model which was assumed to generate chaotic classical dynamics. Here we extend 
and improve these numerical tests by using the conformal billiard model treated above. 

As before we model the quantum dot as a deformed billiard accessible by tunneling from 
leads as shown in Fig. We neglect electron-electron interactions for the following reasons. 
First they will add a charging energy which is irrelevant to the level-width fluctuations. 
Second, although the quasi-bound levels in the presence of electron-electron interactions 
will surely differ from those in its absence, we do not expect this difference to change their 
statistical properties (at least in the chaotic case). This concept underlies the universality 
of RMT and is supported by experimental and theoretical work in nuclear scattering. For 
example complicated shell-model calculations including the residual nuclear interaction lead 
to spectra which exhibit RMT statistic^^. The fact that RMT statistics arise in disordered 
or chaotic non-interacting quantum hamiltonians by no means imply that they are occur 
only when interactions are negligible. 

As in nuclear physicsS one may relate the scattering resonances to the eigenstates 
Xx{u,v) of the dot in isolation using R-matrix theoryll3il3. In the standard approach to 
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elastic scattering from nuclei (for which spherical symmetry may be assumed) a linear rela- 
tionship is derived in each angular momentum channel between the scattering wave function 
and its derivative at the interface between the nucleus and free space. The coefficient of pro- 
portionality is denoted by R{E). If M different decay channels exist this linear relationship 
defines a matrix of coefficients known as the R-matrix for the nuclear reaction. In our case 
we do not have spherical symmetry so an angular momentum expansion is inappropriate, 
but we have a simplification due to the fact that only M propagating modes exist at the 
Fermi energy in each of the two leads. Moreover, in general the tunneling rate (barrier 
penetration factor) will be largest for tunneling into the lowest propagating mode, so we 
can in first approximation neglect all but this mode in each lead. With this approximation 
and taking the simplest case of Neumann boundary conditions for the Xx a derivation very 
similar to that used in the nuclear case yieldslll a 2 x 2 i?-matrix of the form 

l^n{E) = -j:^^, (39) 

A=l ^ - 

where now the index m (= r, /) simply denotes the left or right leads, and the quantities 7]^^ 
and Ex are determined by the solutions of the Schrodinger equation Eq. 

-W/2 
-W/2 



I j'+W/2 

Ix = V^A /_ dn^ y (r/„)XA(C, Vm). (40) 



where Pa is the barrier penetration factor, S,m denotes the direction parallel to lead m and 
r] that perpendicular, Y{rim) is the transverse wave function in the leads (which have width 
W) , and is the position of the inner edge of tunnel barrier m. In deriving this expression 
we assumed that the tunnel barrier is approximately uniform in the transverse direction. If 
we were trying to derive information about specific resonances of this system from the R- 
matrix then we would need to use a more general R-matrix than Eq. (^9]) derived for the GBC 
relevant to our problem. However we are only interested in statistical properties of the R- 
matrix and our results above indicate strongly that the statistical properties are independent 
of the boundary conditions imposed to a good approximation (as long as they allow non- 
vanishing wavefunctions at the tunnel barriers). Therefore we use the more convenient NBC 
for the calculations below. 

An exact non-linear relationship exists between the S-matrix and the R-matrixlil; however 
a particularly useful feature of this formulation is that as E approaches a particular Ex, the 
term in Eq. (|39|) containing the corresponding Xx will dominate (as long as coupling to the 
leads is weak) and all other terms may be neglected. In this approximation the relationship 
between the S and R matrices simplifies to yield the Breit-Wigner formula for the resonance 
line-shape under very general conditions. Moreover the level-width which appears in this 
expression is simply 

= —\iY? ^—p'>r\in\ (41) 

m m 

where we have defined the reduced width 7^. This reduced width will fiuctuate from level 
to level (and for different lead positions for a given level) due to the complicated spatial 
structure of the chaotic eigenfunctions Xx (see Fig. 8). For example, if there happens to be 
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a nodal line near the position of a given lead then the width associated with that level and 
that lead will fluctuate down. (Note that in the absence of spatial symmetry another lead 
attached a few wavelengths away will give a completely uncorrelated width for the same 
level) . 

In order to relate the level-widths to the experimentally-observed peak amplitudes we 
use the expression for the peak height gmax derived by BeenakkerS for the single-level regime 
(assuming kT ^ f ), 

where Tx = F';^ + is the total decay width for level A and F^, F^ are the partial decay 
widths into the right and left leads. The factor ax in Eq. ( ^2]) is a dimensionless measure of 
the area under the T = A-resonance, hence the observed amplitude fluctuations reflect the 
fluctuations in these areas. In this discussion we will only treat the case of equal barriers 
on each side of the dot and hence Px{E) will be the same on the left and right. This means 
that the average decay widths to the left and right are equal, F' = F'" = r/2 and we can 
express ax as 
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The statistics of the peak amplitude fluctuations then follow from those of the reduced 
partial widths |7aP using Eqs. (^) and (|43|) . 

If we assume that the resonance wave functions Xx are described by the GOE when 
time-reversal symmetry is present [B = 0) and by the GUE when time-reversal symmetry 
is broken (sufficiently high magnetic field), then the distribution of partial widths Fx should 
be a xt distribution with the degrees of freedom u = 1 and u = 2 respectivelyli3. This 
distribution should be universal in the chaotic regime, i.e. two different shapes both of 
which generate chaotic classical dynamics should have the same distribution of level widths 
(even though the individual levels are quite different). However if the system approaches 
integrability then non-universal distributions differing from should arise. Precisely this 
behavior is confirmed by our numerical calculations of the partial width distribution for the 



conformal billiard model as seen in Fig. |T2|, where the fully chaotic GOE models (c) and (d) 
fits a xl=i distribution, the fully chaotic GUE models (e) and (f) fits a xl=2 distribution, 
while neither the regular model (a) nor the soft-chaotic model (b) fits a xl distribution. A 
more extensive demonstration of the universality in the chaotic regime is seen in Fig. p!3| , 
where we plot for each value of the shape parameter 6 in the deformation sequence of Fig. ^ 
the value of u yielding the best fit of a xl distribution to the numerically calculated 
histograms both without and with an AB-fiux. It is seen that in all cases the value of u 
fiuctuates around the expected line u = 1 (zero AB-fiux) and u = 2 (non-zero AB-fiux). 

Having confirmed that random-matrix theory works well in the chaotic regime, one can 
derivelll from Eq. (^31) the probability density Vu{<y) where u = 2 for the orthogonal case 
and u = 4 for the unitary case. One finds 



V2{a) = ^2/7ra e"^" (44) 
V^ia) = Aa[Ko{2a) + Ki{2a)] e'^", (45) 
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where Kn are the modified Bessel functions of the second kind. V2 and V4 are plotted in 



Fig. where they are compared to numerical data obtained by evaluating a in Eq. (|43| ) for 
the 300 lowest wave functions of the conformal model for NBC in both the GOE and the GUE 
case. The time-reversal symmetry-breaking needed to study the GUE case is achieved by 
adding an AB-flux of one quarter of a flux quantum. In contrast to the results in Ref. |T2] we 
flnd excellent agreement between random-matrix theory and numerical calculation in both 
the GOE and the GUE case. Note the substantial suppression of small peak amplitudes 
caused by breaking time-reversal symmetry. This reduces substantially the variance of a, 
from Eqs. (|D and (|^) one flnds Aaj/Aal = 32/45 ^ 0.71. 

The effect of a TR-symmetry breaking magnetic fleld on the distribution and its moments 
provides perhaps the simplest experimental test of our theory. However if the suppression of 
the amplitude fluctuations due to time-reversal symmetry breaking is to be cleanly observ- 
able then the magnetic fleld necessary to induce the GOE-GUE transition must be small 
compared to that needed for Landau level formation. Landau level formation strongly sup- 
presses the fluctuation£ll; the classical analogue of this effect is the suppression of chaos 
by the formation of stable skipping orbitJ^i. We estimate the magnetic fleld scale for TR- 
symmetry breaking by adapting an argument flrst put forward (to our knowledge) by Berry 
and RobnikS. TR-symmetry is broken flrst not by the dynamic effect of the fleld but by its 
effect on the phase of the wavefunctions (essentially the Aharonov-Bohm effect). Therefore 
in estimating the TR-symmetry breaking scale we may neglect the dynamic effect of the 
fleld entirely. Gutzwiller's trace formula! implies that structure in the spectrum on the scale 
of the level spacing Ae arises from periodic orbits of period T ^ h/Ae. A magnetic fleld 
will change the action (in units of h) of such orbits by BAT/{h/e), where At is the area 
enclosed by the periodic orbit of period T in the chaotic case. The time-reversed orbit will 
of course enclose area —Ax and their relative phases will be shifted by order unity when 
BAt ~ (h/e), breaking TR-symmetry. Thus the critical fleld Be for TR-symmetry breaking 
is given by Be ~ {h/e)/ At and one need only estimate At- Berry and Robnikil treated the 
case of an A-B flux as above and then evaluated the mean-squared winding number for such 
orbits in the chaotic limit. Their results can be extrapolated to a uniform fleld simply by 
assuming a typical (positive or negative) area of order A (the area of the dot) is enclosed 
with each circuit. With this modest assumption the TR-symmetry breaking flux 

= BeA ~ [Ae^/hvfY'^{h/e). (46) 

Although in the experiments the dot is not isolated as assumed in this argument, the con- 
dition r <^ Ae insures that electrons remain in the dot long enough for the argument to 
still apply. The ratio of this fleld to the fleld at which the cyclotron radius becomes of 
order the radius of the dot scales like N~^^'^, where N is the number of electrons; so in 
dots containing a few hundred electrons TR-symmetry breaking occurs at a fleld one to two 
orders of magnitude smaller than that needed for edge-state formation. In the experimental 
systems of interest this corresponds to a fleld of order a few times 10 mT. Thus the statisti- 
cal effect of time-reversal symmetry breaking predicted by our theory should be observable 
experimentally 
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VII. OPTIMAL EXPERIMENTAL SETUP 



To obtain the most direct experimental verification of tlie tfieoretical results presented 
here and in earlier work it is desirable to fabricate a quantum dot with a variable shape. 
In current experimental systems only a few tens of CB peaks are measured in a given dot 
and these are superimposed on a significant background which complicates the comparison 
of the height of widely-separated peaks. Thus obtaining reasonable statistics (e.g. roughly 
one hundred peak amplitudes) is very difficult. It is possible to use the magnetic field itself 
as an external control parameter causing peak amplitude fiuctuationS, but this only allows 
one to collect statistics for the GUE case. Therefore we suggest an optimal heterostructure 
for tests of RMT statistics would consist of a dot formed by multiple teeth-like gates (see 
Fig. |T3p. By changing the voltages slightly on the inner gates of this structure it should be 
possible to realize a whole range of shapes without affecting the region near the quantum 
point contacts defining the tunnel barriers to the surrounding two-dimensional electon gas. 
It would then be possible to follow a given peak as a function of shape and collect statistics 
without the complication of background variation, both in the presence and absence of a 
magnetic field. Since the magnetic field would not be the control parameter it would also 
be possible to map out the GOE to GUE transition. This concept is illustrated by the 
calculations for the deformed billiard shown in Fig. 16 where we plot the peak amplitude 
a for two given levels as function of the shape parameter b. A particularly simple quantity 
to extract from such experiments is the correlation function for the peak amplitude as the 
shape is varied. Although we are not able to make quantitative estimates at this point our 
results suggest that for dots containing a few hundred electrons (as in current experiments 
on micron-size devices) quite small voltages would be needed to rearrange the wavefunctions 
of the highest levels and hence decorrelate the peak amplitudes. In the single-level regime 
[kT < Ae) the theory predicts that the voltage scale should be independent of temperature; 
a prediction which could be experimentally tested and (if confirmed) would provide support 
for the basic model. 



VIII. DISCUSSION AND CONCLUSION 

In this paper we have studied the statistical properties of deformable billiards in the 
mixed and fully chaotic limits, with applications suggested to quantum dot systems. We have 
presented numerical results based on the continuous family of conformal quantum billiards 
introduced by Robnik for which efficient numerical solution of the Schrodinger equation is 
possible and for which a rather complete characterization of the classical mechanics exists, 
including the existence of a KAM transition. 

This model system has provided us with a unifying framework within which we discussed 
several aspects of quantum chaos in billiards. Our calculations demonstrate the universal 
behavior of the system once chaos is sufficiently strong (hard chaos) in which case random- 
matrix theory describes the quantum statistical properties independent of the details of the 
billiard geometry or the boundary conditions on the quantum problem. The calculations 
also demonstrate that in the mixed regime (soft chaos) the behavior is not universal and the 
random-matrix ensembles do not describe the system. We have also verified the theoreti- 
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cal results obtained in the literature for the parametric energy level correlations of isolated 
quantum dots and the statistics of Coulomb blockade conductance peaks based on the as- 
sumption that RMT or the supersymmetric a model applies to chaotic billiards thereby 
adding credibility to these assumptions. Finally, we have proposed an experimental setup 
which is particularly well suited for tests of the theoretical predictions in this work. 
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FIGURES 



FIG. 1. The deformed billiard T) is the image of the unit disk C under the mapping w{z), 
which is conformal {w'{z) 7^ 0, Vz G C) and area preserving (/^ \w'{z)\^dz = vr). 

FIG. 2. Two sequences, (a) and (b), of dV as a function of changing parameters. Sequence 
(a) was studied by Robnik0. Sequence (b) is studied in this paper. For the very last shape in each 
sequence w{z) is in fact not conformal. In both of these extreme cases w'{z) = just at the point 
z = -1. 

FIG. 3. Four Poincare surface of sections for the deformation sequence shown in Fig. 

FIG. 4. The Kolmogorov entropy plotted versus the normalized deformation parameter X for 
the two deformation sequences (a) and (b) in Fig. ^. In (a) X = b, and in (b) X = 5. 

FIG. 5. The level spacing histogram for the Africa billiard b = c = 0.2 and 6 = tt/S with an 
AB-flux of and with DBC. A x^-test supports the hypothesis of GUE level statistics (the full 
line) (P(x^) = 0.68), whereas the hypothesis of GOE level statistics (the dotted line) is rejected 
(P(x2) = 0.0002). 

FIG. 6. The t^st for the deformation sequence shown in Fig. |2|b in the case of DBC and 
with an AB-flux of j^o- The full line (*) is the test for GUE statistics while the dashed line (o) 
is for GOE statistics. The horizontal line is the (logarithm of the) mean value 1/2 around which 
a successful -test ought to fluctuate. Note the cross-over between the two statistics near the 
symmetry points 6 = and tt. 

FIG. 7. The test for the deformation sequence shown in Fig. |2|b in the case of NBC and 
with an AB-flux of j^o- The full line (*) is the test for GUE statistics while the dashed line (o) 
is for GOE statistics. The horizontal line is the (logarithm of the) mean value 1/2 around which 
a successful -test ought to fluctuate. Note the rather wide cross over region between the two 
statistics near the symmetry point 6 = 0. 

FIG. 8. A deformed billiard with two leads attached. The leads can either be completely open 
or they can connect to the billiard through a tunnel barrier. The natural boundary conditions in 
this case are GBC. The wavy lines inside the billiard is the nodal structure of an eigenstate as 
briefly discussed in Sec. VI. 

FIG. 9. The x^ t^st for the deformation sequence shown in Fig. |2|b in the case of GBC and 
zero AB-flux. The full line (*) is the test for GUE statistics while the dashed line (o) is for GOE 
statistics. The horizontal line is the (logarithm of the) mean value 1/2 around which a successful 

test ought to fluctuate. 
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FIG. 10. Part of the normalized energy spectrum (levels 225 to 274) of the deformation 
sequence Fig. |2|b with an AB-flux of one quarter of a flux quantum and with DBC. 



FIG. 11. In (a) is shown the numerical calculated energy level correlation function c(x) for 
the fully chaotic deformation sequence Fig. ||b both with (*) and without (o) an AB-flux. The 
full (dotted) curve is the universal correlation function calculated with (with out) an AB-flux for 
disordered systems in Ref. 24. In (b) c{x) is shown with zero AB-flux in both the soft chaos regime 
(A) of the first half of the deformation sequence Fig. |2|a and the hard chaos regime (o) repeated 
from (a). Note the strong dip for small values of x in the case of soft chaos. 



FIG. 12. The distribution of the partial decay width, Fx, for six different deformed quantum 
billiards. Shown are histograms representing our numerical results, xt distributions (full curves) 
serving as a guide for the eye = 1.0 in (a)-(d) and 2.0 in (e)-(f)), and as inserts the particular 
shapes of the dot. In (a), (b), and (c) we used the simple quadratic map with c = and with b = 
0.00, 0.14, and 0.39 respectively. In (d), (e), and (f) we used the cubic map with b = c = 0.2 and 
6 = vr/3, 27r/3, and tt/S . In (a) the model is integrable, in (b) it is nearly integrable, and in (c)-(f) 
it is fully chaotic. There is no AB-flux in (a)-(d) and one quarter of a flux quantum in (e) and (f). 



FIG. 13. The values of obtained from fits of a Xu distribution to the numerically calculated 
histograms of the distribution of the partial decay widths for the entire deformation sequence 
in Fig. Qb with zero AB-flux (o) and with an AB-flux of l/4<I>o (*)• 



FIG. 14. Predicted distribution of peak amplitudes a in the presence (a) and in the absence (b) 
of time-reversal symmetry, compared to the numerically generated amplitude distribution obtained 
with the shape parameters 6 = c = 0.2 and 6 = tt/3. 



FIG. 15. A top view of a gated heterostructure. The gate is splitted into twelwe teeths. The 
two pairs at the ends form quantum point contacts leading to the surrounding two-dimensional 
electron gas. The remaining eight teeths define the shape of the quantum dot. In situation (a) 
all the interior gates have the same voltage. In (b) one gate has a slightly higher voltage and 
another slightly lower thereby deforming the dot while maintaining its area. The regions close to 
the quantum point contacts are essentially unaffected by this action. 



FIG. 16. The variation in the peak amplitude a for two particular levels as a function of shape 
deformation 6. 
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